arXiv:1507.01940vl [astro-ph.CO] 7 Jul 2015 


Self-similarity of Dynamo Action in the Largest Cosmic 
Structures 


Francesco Miniati^ & Andrey Beresnyak^ 

^Physics Dept, ETHZurich, Wolfgang-Pauli-Strasse 27, CH-8093 Switzerland 

‘^Nordita, KTH Royal Institute of Technology and Stockholm University, SE-10691 Stockholm, 

Sweden 

Massive galaxy clusters (GC) are filled with a hot, turbulent and magnetised intra-cluster 
medium (ICM). Still forming under the action of gravitational instability they grow in mass 
by accretion of supersonic flows. These flows partially dissipate into heat through a complex 
network of large scale shockiP, while residual transonic flows create giant turbulent eddies 
and cascadd^l^ Turbulence heats the ICiyP and also amplifies magnetic energy by way of 
dynamo actionP®. However, fundamental properties of the pattern whereby gravitational 
energy turns kinetic, thermal, turbulent and magnetic remain unknown. Here we find that 
the energy components of the ICM are ordered according to a permanent hierarchy, in which 
the ratios of thermal to turbulent to magnetic energy densities remain virtually unaltered 
throughout the ICM history despite evolution of each individual component and the drive 
towards equipartition of turbulent dynamo. Our results are based on a state-of-the-art, fully 
cosmological computational model of ICM turbulenc€013, revealing that an approximately 
constant efficiency of turbulence generation from gravitational energy that is freed during 
mass accretion. The permanent character of this hierarchy reflects a new aspect of self- 


1 


similarity in cosmologj®®', while its structure, consistent with current data ISHSl, encodes 
information about the efficiency of turbulent heating and dynamo action. 

The computational model captures the turbulent motions through a multi scale technique 
which employs six nested grids covering progressively larger volumes with correspondingly coarser 
resolution elements The finest grid resolves the virial volume of the GC with more than a bil¬ 
lion uniform-size resolution elements and provide the necessary dynamic range to resolve the ICM 
turbulent cascade. The largest grid covers the chosen cosmological volume of 340 comoving Mpc 
on a side (comoving=partaking in the expansion of the universe; 1 Mpc ~ 3 million ly). The 
intermediate grids allow to simultaneously follow with adequate accuracy the matter distribution 
outside the GC volume. The calculation starts with three grids and adds progressively finer grids 
as the Lagrangian volume of the GC shrinks under self-gravity. All six grids are in place at a time 
corresponding to 8 billion yr after Big Bang. At current time (13.8 billion years after Big Bang) 
the simulated GC has a total virial mass of 1.3 x 

Figure shows a snapshot of the simulation illustrating the cosmological context and the 
highly turbulent conditions of the flow inside the GC volume. The exquisite resolution across the 
GC volume allows us to accurately measure the time dependent statistical properties of structure 
formation driven ICM turbulence including, in particular, the dissipation rate, eturh, the outer scale, 
L and the velocity dispersion on that scale, (Methods). In the following we restrict 

our analysis to a region within 1/3 of the GC’s virial radius, Rvir, where Rmr defines a region 
characterised by a mass over-density Ac ~ 100 that has nominally reached dynamical equilibrium. 
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Our choice is motivated by the faet that at eurrent epoch Rmr/3 ^ 1 Mpc, i.e. defines a region 
most relevant for eomparison with observations. 

Hydrodynamie turbulenee is dominated by the solenoidal eomponent aeeounting for 60-90% 
of the total kinetie energjEE^. Detailed analysis shows that this eomponent remains statistieally 
homogenous and isotropie thus resembling Kolmogorov’s caseade, despite the presenee of eon- 
siderable strueture in the IClvP. The dissipation of ineompressible turbulenee eontributes to ICM 
heating along with shocks and adiabatic compression and to the growth of magnetie energy by 
way of small seale dynamo actiorP^ (Methods: Fig. |^for caseade details). The turbulent dis¬ 
sipation rate associated to the solenoidal eomponent is estimated from the numerieal simulation 
data. Beeause ICM turbulenee is driven by various complex hydrodynamie meehanisms ultimately 
powered by the unsteady mass aeeretion proeesP, the dissipation rate is highly ehangeable with 
time and exhibits non-monotonie variations by more than one order of magnitudeP' (Figure]^). 
However, alongside with mueh eomplexity turbulent dissipations appears to also exhibit simplieity 
of behaviour. This is shown in Figure illustrating the time evolution of the fraction of thermal 
energy originating from turbulent dissipation. In contrast to eturb, this quantity remains remarkably 
eonstant during the GC lifetime, rjturb ~ 0.3 — 0.4, indieating that the effieiency of turbulenee gen¬ 
eration out of gravitational energy freed by mass aeeretion is approximately eonstant. In addition. 
Figure shows that the turbulence velocity dispersion at the outer seale normalised to the ICM 
sound speed, i.e. the turbulenee Mach number M.turb = ^ /Cs, also remains rather eonstant 

with time. This shows that in the ICM the evolution of the turbulent kinetic energy and the thermal 
energy are elosely related, eonsistent with the previous plot. The value of Mturb can be under- 
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stood as follows. If the generation of the bulk of the thermal energy, Et^, is dominated by the last 
a = 2 - 3 eddy turnover times, tl = L/ul, then Eth ^ peturbdt ~ {a/3^turb)p{{5uLf), 
where we have used the known relation eturb = (2/3(7) i {{SulY) i/L with C ^ It is straight¬ 
forward to then see that the Maeh number ~ {V3/a)Yvturb/0.37)^ whieh,forQ! = 1.5 — 3 
ranges between 0.8 and 1.2 eonfirming the result in Figure It also follows that the ratio of 
thermal to turbulent kinetie energy is 

dd'th 2tt 

IpHSulY) si 

Generation of magnetie field by small seale dynamo in a turbulent flow follows from standard 
theory. In a high Reynolds number (Re > 10^) flow sueh as the ICM an initial seed of vanishing 
strengthP^^^ is amplified exponentially at the rate 7 = y/^/SOri, where = L/5ul is the eddy 
turnover time. After a short while (oc magnetie field stops growing below a eharaeteristie 

Alfven seale. La = v\IC^eturb, where va = S/^/dvrp is the Alfven speed, due to the feedbaek 
aetion of magnetie tensiorl^Il. Magnetie energy eontinues to grow at the expenses of turbulent 
kinetie energy as La, marking the equipartition seale between kinetie and magnetie energy, shifts 
towards larger valued. Growth, however, is now proportional to the turbulent dissipation rate 
instead of exponential with time. It is in this latter stage that the dynamo spends most of the 
timeP. Reeent state-of-the-art numerieal work finds that for statistieally isotropie and homogeneous 
turbulenee, as found in the IClVi 2122 l, the effleieney of eonversion of turbulent (kinetie) to magnetie 
energy is a universal number ea Ce= 4-5%P. 

Therefore, the evolution of magnetie energy in the ICM ean be expressed in terms of the tur- 
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bulence dissipation history as -Es(t) = = Ce /* dTpe{T). Combined with the above finding 

about rjturb, this leads to simple but significant expressions relating the fundamental properties of 
magnetic field and and turbulence in the ICM. In fact, since turbulence dissipation contributes a 
constant fraction, pturb ~ 1/3, of ICM thermal energy, the ratio jdpiasma of thermal pressure to 
magnetic energy can be written as 
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This means that for massive GC fdpiasma is a constant, which depends neither on the specifics of the 
ICM conditions including turbulence, nor on the GC mass or age. It is instead simply determined 
by two fundamental parameters, Ce and rjturb, which describe the efficiency of turbulent dynamo 
and of turbulent heating in structure formation, respectively. This is shown in Figure where 
is plotted as a function of cosmic time and exhibits 25% rms fluctuations, which should 


also characterise massive cluster-to-cluster variations. We can also compute the Alfven scale. 
Since the turbulence is non-stationary and the magnetic energy retains memory over more than 
one eddy turnover time, we average the dissipation rate eturb over 2 Gyr when calculating La- 
Expressing La in units of the turbulence outer scale we write 


La ^ v\ ^ 3 / 2 \ ^ 

L C"^ {eturb) 2 \^(dplasma) {{Sul)"^)^ 


1 ^ ^ 

100 


40 




( 3 ) 


We have already shown that both /3piasma and the turbulence Mach number, Mturb = ^/ Cs, 

remain constant during the evolution of the GC. Therefore, the Alfven scale too remains a constant 
fraction of the turbulence driving scale, independent of time, GC mass and ICM conditions. In 
addition, given the large value of /dpiasma, and that Mturb ~ 1, is small compared to L. The 
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time evolution of is shown in Figure]^ (see Fig. |^in Methods for a typieal ICM spee- 
trum of hydromagnetie turbulenee). Finally, Figure]^ shows that the evolution of the turbulent 
injeetion seale, L, elosely follows that of La while traeking the growing eharaeteristie seale of 
the GC (i? 5 oo = also plotted in the same panel. Note that the modulation of La refleets 

the ehanging turbulent eonditions in the ICM and, in partieular, is anti eorrelated with eturb, as 
generally expeeted. 

The large value of (3piasma indieates that magnetie energy, like turbulent energy, is small 
eompared to thermal energy {Eth S> Eb). Moreover, the small value of La/L indieates that the 
dynamo is far from saturation and magnetie energy is also small in eomparison to the turbulent 
kinetie energy {Eturb 3> Eb). This energy hierarehy is fundamentally due to the effieieney T]turb 
with whieh turbulent energy is generated during gravitational eollapse and the fraetion Ce thereof 
that is eonverted into magnetie energy, namely Eth ■ Eturb ■ Eb = I '■ rjturb ■ CEVturb- The values 
of /Spiasma and La are in good agreement with reeent measurements of magnetie field properties 
in GC nSUll Here, they emerge from pure numerieal modelling of strueture formation turbulenee 
and MHD dynamo aetion, in the sense that the are found to derive their values from the parameters 
r]turb and Ce, whieh are determined numerieally and not through parametrie fits. Intriguingly, the 
above energy hierarehy appears to remain unehanged during the GC evolution and the turbulent 
dynamo in the ICM is far from saturation today as it has virtually always been in the past. Figure]^ 
shows that the GC size and, therefore, its mass eonstantly grow. This implies that the gravitational 
potential energy and therefore the ICM thermal energy and turbulent energy also eontinue to grow. 
Meanwhile dynamo aetion tries to bring magnetie and turbulent energy into equipartition. Sinee 
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all of these forms of energy grow simultaneously but with different eonstant effieiencies, their 
ratio remains unehanged, reflecting the value of those intrinsic efficiencies. In other words, both 
Ppiasma ^nd Lyi/L cncodc the efficiency of turbulent generation in structure formation and the 
efficiency of dynamo action. As such, they allow us to relate magnetic field observations in massive 
galaxy clusters to such properties of structure formation. This is in sharp contrast with other 
astrophysical bodie^^^®^, e.g. the interstellar medium of galaxies, stars and compact objects, where 
the turbulence dynamo has long saturated and such information is lost forever. 
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Methods 


Numerical model The simulation is carried out with CHARM, an Adaptive-Mesh-Refinement cos¬ 
mological codeia. This code uses a directionally un-split variant of the piecewise parabolic method 
for hydrodynamic^, constrained-transport algorithm for solenoidal MHDI^SI^ a time-centred mod¬ 
ified symplectic scheme for the collision-less dark matter, and solve Poissons equation with a 
second-order accurate discretisation. The magnetic field remains negligible throughout, so the 
calculation is effectively hydrodynamic. For massive galaxy clusters, such as Coma cluster, the 
ICM cooling time is a few times the age of the universe!^, so cooling and baryonic feedback 
processes are neglected. Heating of the intergalactic medium through photoionization is also ne¬ 
glected, with no consequences whatsoever for the generation of vorticity and turbulence at ac¬ 
cretion shocks. We use a concordance A-CDM universe with normalized (in units of the critical 
value) total mass density, = 0.2792, baryonic mass density, Vti, = 0.0462, vacuum energy den¬ 
sity, = 1 — Vtm = 0.7208, normalized Hubble constant h = Hq/IOO km s“^ Mpc“^ = 0.701, 
spectral index of primordial perturbation, Ug = 0.96, and rms linear density fluctuation within a 
sphere with a comoving radius of 8 h~^ Mpc, as = 0.817®. The simulated volume has comoving 
size of Lbox = 240 h~^ Mpc on a side. The initial conditions are generated on three refinement 
levels with graf ic++ (made publicly available by D. Potter). For the coarsest level we use 512^ 
comoving cells, corresponding to a nominal spatial resolution of 468.75 h~^ comoving kpc and 
512^ particles of mass 6.7 x 10® h~^ Mq to represent the collisionless dark matter component. The 
additional levels allow for refined initial conditions in the volume where the galaxy cluster forms. 
The refinement ratio for both levels is, nfgf = = 2, f = 0,1. Each refined level covers 
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1/8 of the volume of the next eoarser level with a uniform grid of 512^ eomoving eells while the 
dark matter is represented with 512^ partieles. At the finest level the spatial resolution is Ax = 
117.2 h~^ eomoving kpe and the partiele mass is 10® h~^ Mq. As the Lagrangian volume of the 
galaxy eluster shrinks under self-gravity, three additional uniform grids eovering 1/8 of the volume 
of the next eoarser level are employed with 512®, 1,024® and 1,024® eomoving eells, respeetively, 
and nfgf = 2,4, 2, for f = 2,3,4, respeetively. All of them are in plaee by redshift 1.4, providing a 
spatial resolution of 7.3 h“® eomoving kpe in a region of 7.5 h~^ Mpe, aoeommodating the whole 
virial volume of the GC. The ensuing dynamie range of resolved spatial seales is suffieiently large 
for the emergenee of turbulenee. 

Galaxy cluster characteristic quantities The galaxy eluster and its formation history are reeon- 
strueted using our implementation of a HOP halo findei^ and merger history eode. The virial 
radius is defined as the region enelosing a mass over-density Ac = with respeet to the 

eritieal densit}^^. At redshift 2 ; = 0, the viral radius is R^r = 1.95/i“® Mpe, and the eorresponding 
enelosed mass, M^ir = 1.27 x Also a z = 0, using Ac = 500 we find the eharaeteristie 

radius i /500 — Mpe. 


Turbulence characterization The eharaeteristie quantities deseribing the turbulenee are inferred 
from the analysis of the strueture funetions. This analysis is deseribed in detail irl^H. Basieally, we 
deeompose the veloeity into a solenoidal and a eompressional eomponent using a Hodge-Helmoltz 
deeomposition, i.e. 


V = -f Vc, Vc = - V0, = V X A, 0 
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and then we eompute the seeond and third order strueture funetions of veloeity inerements of the 
solenoidal eomponent, 5vi = [vs(x + 1) — Vs(x)]i^, 

^,(1) ^ , (5) 

where p = 2, 3 indieates the strueture funetion order, and i indieates the projeetion along or per- 
pendieular to 1 for the longitudinal and transverse strueture funetions respeetively. To eompute 
the strueture funetions we define sampling points randomly distributed inside the volume of in¬ 
terest (within (1/3) of the viral radius), and eompute the veloeity differenee with respeet to other 
randomly seleeted field points at a maximum distanee of two virial radii. Onee we the veloeity 
strueture funetions are eomputed, we define the veloeity dispersion as the asymptotie values of 
the seeond order strueture funetion, and the outer seale as the separation at whieh that asymptotie 
value is reaehed. To eompute the Maeh number we divide the turbulent veloeity dispersions by 
the sound speed, Cg = \JlPlp, eomputed by evaluating the mean value of eaeh thermodynamie 
quantity within the same volume in whieh the sampling points are eolleeted. Finally, the turbulent 
dissipation rate is eomputed by identifying the inertial range of the seeond and third order strueture 
funetions of the solenoidal veloeity inerements (for details see ref.^^. 

Code Availability We have opted not to make the eode available for praetieal reasons. However, 
the methods we adopt are published in the literature and are eommonly used in the eommunity. 
Amongst others, the publieations mentioned in the above Methods seetion eontain tests of our eode 
against problems with known solutions and also with respeet to solutions obtained with similar 
eodes from independent authors. 
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Figure 1 | High resolution simulation of a galaxy cluster in fully cosmological context. Bary- 
onic gas in the large scale structure of the universe (panel a; bright is high, dark is low), 
and around the GC centre where numerical resolution is highest (panel b; colormap in¬ 
verted). Dynamic range of density (in cm~^) is ^ 10®. The black dash-line marks the 
virial radius, R^ir, enclosing the volume that has nominally reached dynamical equilib¬ 
rium. Panel c: vorticity magnitude on a scale twice the finest mesh size (~ 20 kpc). 
Panel d: temperature map. Complexity is due to shocks and contact discontinuities in the 
turbulent flow. Dynamic range of temperature (in K) and vorticity (in is about 10^. 

Figure 2 | Time evolution of turbulence. Panel a: rjturb, the ratio of ICM thermal energy 
contributed by turbulent dissipation, J^dtpeturb, to the total thermal energy, Eth- Panel 
b: turbulent Mach number, the ratio of the turbulent rms velocity to the sound 

speed Cs. Panel c: volumetric turbulent dissipation rate, eturb obtained irP. The error bar 
correspond to the variance of eturb- All quantities are computed within 1/3 of the virial 
radius. Time in billion-yr is reported on the bottom x-axis and cosmological redshift on the 
top x-axis. 

Figure 3 | Time evolution of magnetic field. Panel a: Ppiasma, the ratio of ICM thermal to 
magnetic pressure computed as B'^/Sn = Ce f dTpeturb{r). Panel b: La/T, the ratio of 
Alfven to the turbulent injection scale. La{t) = v\/[C^eturb)], where va = B/y/^ is the 
Alfven speed and (eturb) is the turbulent dissipation rate smoothed over r = 2 Gyr with 
a Gaussian filter. Quantities in panels a,b refer to a volume inside 1/3 the virial radius. 
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The dash lines show transients to the asymptotic regime for an artificial tstart = 4.5 Gyr. 
Panel c: turbulence injection scale L (solid-line) and the characteristic cluster size R^oo = 
0.5 R^ir (dash-line), enclosing a mass over-density of 500. Time in billion-yr (bottom x- 
axis) and cosmological redshift (top x-axis) are reported. 

Figure 4 | Generation and cascade of ICM hydromagnetic turbulence. First the gravita¬ 
tional potential energy is converted into kinetic energy of accretion flows. These generate 
shear and shocks which, in addition to heat dissipation, produce fluid instabilities and 
baroclynic term, respectively, leading to turbulent flows. Shocks also accelerate particles 
through the Fermi I mechanism. Shocks do not dissipate tangential flows which will either 
generate turbulence, shear or shocks or a combination thereof. The turbulence cascade 
includes, dissipation of compressible modes at weak shocks, conversion of turbulent to 
magnetic energy via dynamo action, excitation of plasma waves accelerating relativistic 
particles through Fermi II mechanism, and of course viscous dissipation. 

Figure 5 | Spectrum of ICM hydormagnetic turbulent cascade. Characteristic spectrum of 
turbulent kinetic energy in the ICM. Solid and dashed lines correspond to the solenoidal 
(Kolmogorov-like) and the compressional (Burgers-like) velocity field, respectively. On the 
X-axis, from left to right we have marked the virial scale, R^ir, the injection scale, L, the 
Ozmidov’s scale, Lq, the Alfvn scale. La, and Kolmogorov’s dissipation scale, idiss- All 
quantities are time dependent and Ozmidov’s scale is comparable to the injection scale, 
so at times turbulence in the radial direction could be suppressed by stratification. 
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